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Abstract. Recent advances in the treatment of scattering of charged compos- 
ite particles are reviewed. In a first part I report on developments of the theory. 
Specifically I describe the recent completion of the derivation of the co-ordinate 
space asymptotic behaviour of the wave function for three charged particles in 
the continuum. This knowledge is increasingly being made use of in attempts to 
'derive' three-Coulomb particle wave functions to be used in all of configuration 
space which are solutions of the Schrodinger equation, though not everywhere 
but at least in one or preferably all of the asymptotic regions. Their practical 
application in approximate calculations of ionisation and breakup processes is 
pointed out. The asymptotic three-charged particle wave functions find further 
use in investigations of asymptotic and analytic properties of matrix elements 
of the three-body Coulomb resolvent. An important example is the nonper- 
turbative derivation, valid for all energies, of the large-distance behaviour of 
the optical potential. In the second part I describe a renewed attempt to es- 
tablish the few-body approach as a valuable tool for calculating (energetic) 
atomic collision processes. At the end I briefly touch upon the recent successful 
above-breakup-threshold calculation of proton-deuteron elastic scattering for a 
realistic potential. 



1 Theoretical Developments 

1.1 Asymptotic Wave Function for Three Charged Particles in the Continuum 
1.1.1 Introduction 

Knowledge of the asymptotic form of the wave function is required, e.g., as 
boundary condition for selecting the physical solution of the Schrodinger equa- 
tion (SE) for energies above the ionisation threshold, and for deriving asymp- 
totic and/or analytic properties of quantities which are matrix elements of the 
three-body Coulomb resolvent, such as the optical potential or the kernel of 
momentum space integral equations for charged particles (for a discussion of 
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the latter topic, see [jj]). And, last but not least, it avails also for approximate 
calculations of breakup or ionisation in peripheral collisions. 

I start by introducing some notation. Consider three distinguishable par- 
ticles with masses m u and charges e„, v = 7. In the cm. system Ja- 
cobi co-ordinates are used: k a (r a ) is the relative momentum (co-ordinate) 
between particles (3 and 7, and q a (p a ) the momentum (co-ordinate) of parti- 
cle a relative to the center of mass of the pair 7). Frequently, two Jacobi 
three-vectors are combined into a sixdimensional vector: P := (fc„, q v ) and 
X := (r„, p v ), v = a,/3,j. 

For three particles in the continuum the energy shell condition reads E = 
q^/2M a + k^/2p a , with the standard definitions p~ x = m^ 1 +m~ 1 and M" 1 — 
771 q 1 + (mp+'m~ l )~ 1 of the appropriate reduced masses. The interaction between 
particles (3 and 7 is assumed to consist of a short-range and a Coulomb part 

V a = V% + VZ, with V C(r a ) = ^. (1) 

For the following it proves convenient to eliminate the three-body plane from 
the wave function by introducing a 'reduced' wave function #j, +) (X) via 

V [ +\X) =: e ipx ¥+\X) ee e i(k«^+ 9 „^)^ (ra)Pa ), (2 ) 

1.1.2 Asymptotic Solution of the SE in fl®: 'Redmond Asymptotics'. 

This region of configuration space is reached by letting r a , rp , and r 7 ap- 
proach infinity, but such that not r v j ' p v — > for = a, (3, or 7; in other 
words, asymptotically all inter-particle distances grow uniformly. Introduce the 
Coulomb parameter, e.g., for the subsystem 7) as r/ Q = efje 1 p a /k a . Then, 
as stated by Redmond [|| , in leading order one has 

3 

i/=l 

in the so-called 'non-singular regions' (i.e., for k v r v ^ k„ ■ r u ). An asymptoti- 
cally equivalent form is 

*£\X) ~° n^(r„), (4) 
i/=i 

where ipk a (fa), when multiplied with the appropriate plane wave, is the con- 
tinuum solution of two-body SE 

{It + |S " ¥ « {ra) ~ v « {ra) ] ^'^(r*) = °> ( 5 ) 

with leading asymptotic behaviour 

^fe +) (^a) r ^° e i " aln ^ r <"- fe "- T '«). (6) 



For purely Coulombic interactions this is just the known asymptotic form of 
the explicit solution of (||), 

^{r a ) = ^l + \r a )=N a F{-\7 la ,l ] i{k a r a -k a -r a )) (V^ = 0). (7) 

Here, N a = e'^ 2 T{1 +iry a ); F(a, b; x) is the confluent hypergeometric, and 
r{z) the Gamma function. In the general case the behaviour (j|) follows from 
the fact that for sufficiently large r a , v£ (r a ) can be neglected in (|1) yielding 
r ^—>° (T a ). I, finally, mention that the form (Q) with purely 

Coulombic wave functions (^) has been proposed in ||, ||. 

1.1.3 Asymptotic Solution of SE in v — a, (3 or 7: 'AM Asymptotics'. 

Besides J?o, there exist three more asymptotic regions in which the boundary 
condition on the wave function must be prescribed. They are characterised by 

■ Pv — > 00, r u /p u — > 0, for v = a,(3, or 7. (8) 

That is, e.g. in J? Q , the distance p a of particle a from the center of mass of the 
pair ((3, 7) grows faster than the relative separation r a between the particles of 
this pair. Recall the asymptotic behaviour (^) in i?o which we rewrite as 

This suggests that to go from J?o to Q a , which includes all values r a £ [0, 00) 
while still rp, r 7 — > 00, it might suffice to replace in (^) the asymptotic wave 
function (||) for the pair ((3,j) by the full solution $fe a (r a ) of the SE (||). 

However, this simple idea turns out to be not correct. Instead, as has been 
shown by Alt and Mukhamedzhanov Q (AM) the leading term is given by 

& €(pj^) n^ 111 ^^ 7 ^ ( io ) 

(for further and next-to-leading order contributions see || 0, ||, [|). What 
appears is the exact continuum solution of two-body- like SE (for all r a £ fi a ) 

^ + " K?(ra) - Kf (r a )} e*-^-- 4» o) (r o) = 0. (11) 

It describes the relative motion of p and 7 at an energy which depends paramet- 
rically on the distance from particle a, as determined by the local momentum 

MpJ=fca + ^, a a (p a )=- J2 VuKv' e ™P<*- k » , (12) 

Here, A„„/ = m„/(m„ + m„<), fc„ := k„/k v , e a p = —cp a = +1 if (a,/3) is a 
cyclic ordering of (1,2,3). If = 0, the explicit solution of (11) is 

^fel(pj( r «) = ^(Pa) ^(-"^(pJ^KMPaK - k a { Pa )-r a ), (13) 
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which is of the same form as (Q) but with the asymptotic momentum k a re- 
placed everywhere by the local momentum k a (p a ). 

The physically decisive point is that the parametric dependence of the 
wave function V^ + j p ) f° r the pair 7) on the relative co-ordinate p a be- 
tween its center of mass and particle a is a manifestation of a long-ranged, 
noncentral, velocity-dependent, dynamic three-body correlation which is a new 
genuine three-body effect. This is evident from rewriting the bracket in as 



k- 2 
2p a 



2p, a 



V a (r a ) 



1 a a (p a )-k a 



//-n 



o 



Pa 



(14) 



Comparison with the bracket in the ordinary two-body SE (|^) reveals that 
the additional, potential-like term in (14) is indeed of long range (~ Pq 1 ), 
and depends on the relative velocity k a /p a within the pair (/3, 7) and on the 
orientation of k a relative to p a (cf. the definition of a a in jl^)). The latter 
properties can intuitively be understood as a result of the 'dynamic screening' 
of the charges when the three particles are in typical configurations pertaining 
to the asymptotic regions Q v , v — a,/3, / y. 

When can one expect to observe effects of these three-body correlations? As 
follows from (|lj), local (k a (p a )) and asymptotic momentum (k a ) will differ 
appreciably only if k a is small (recall that p a — > 00 in J? Q ). As a measure of 
the 'size' of their difference we can take |a Q (p Q )| which for k a — > behaves as 



Ma 



(15) 



That is, it is inversely proportional to the relative velocity of particle a and 
the center of mass of the pair (/3,7), and proportional to the famous charge- 
over-mass ratio difference of the pair (f3,j) ('Post-acceleration effect'). 



1.1.4 Asymptotically Correct Three-body Coulomb Wave Functions 

Recently many attempts at 'deriving' three-body Coulomb wave functions, 
which extrapolate the known asymptotic wave functions into the non-asymptotic 
region, have appeared. Obviously, such are highly non-unique and, thus, open 
to further optimisation. For instance, a wave function which is asymptotically 
correct in all regions f2\ U f2 2 U i7 3 U i? has been proposed in H 

3 

S { + } {X) = 11 N v F(-ifj„, 1; i(k v r v -k„-r v )), (16) 

with k a = k a + 2a a (p a )/R, R = J2l=i r vi Vf : = Vu(k v ), N u := N v (fj v ). A 
structurally identical wave function but with a different definition of the local 
momenta has recently been proposed in [ pToj, p!T[ . For a wave function which 
contains - again differently defined - local momenta in only two hypergeometric 
functions and hence is asymptotically correct in three regions, say fl a U ftp U , 
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only see |Q . It is to be noted that all these admissible choices of local momenta 
differ from the original one fa] in next-to-leading order in 1/ p a only. 

Most effort has, however, been directed towards devising wave functions 
which are asymptotically correct in J?o only. The prototype is (Q) (with — > 

V'fe ^) w hich, though it has the unphysical feature that the motion of the three 
particle pairs is completely uncorrelated in all space, is nevertheless widely 
used. To overcome this deficiency correlations are introduced, by modifying the 
Coulomb parameter which occur in each of the hypergeometric functions (R): 
r/ a — » /i Q [e^e 7 + x Q ,]/fc Q where either \a — Xa(k a , q a ) depends on the momenta 
('velocity-dependent charges' jO, or on the positions Xa = Xa(r a ,P a ) 
('position-dependent charges' 15 ). The fairly arbitrary functions Xa are chosen 
such that the resulting three-body wave functions satisfy additional criteria like 
the cusp condition. For very general wave functions of this type see [O . 



1.1.5 Applications 

These - at least asymptotically correct - wave functions are increasingly being 
applied to approximate calculations of ionisation and breakup processes. In its 
simplest form (with incoming plane wave \q a )) the corresponding amplitude is 

T amfi (q a ;k' ,q' ) = {^] q ^\V a \i, am )\q a ) , (17) 

with V a = ^2 v ^ a V v being the (entrance) channel interaction and ip a m = 
!/i Qm (r tt ) denoting the incoming bound state wave function. Expression (O) 
contains for the final state the full solution (<?v, , I of the three-charged parti- 

K l3'1f3 

cle SE which is being approximated by one of the aforementioned ansaetze. 

Many calculations have appeared for (e,2e) reactions on hydrogen and other 
atoms; for a comparative study of various types of asymptotically correct wave 
functions see [^6|. A nuclear physics example is the Coulomb dissociation of a 
projectile in the field of a fully stripped nucleus which is an important tool for 
obtaining astrophysical S- factors p7[ . Quite generally one expects the better 
results with such an approximation the more peripheral the collision is (the 
impact parameter should be larger than the range of the nuclear interaction). 



1.2 Long-range Behaviour of the Optical Potential (OP) 
1.2.1 Introduction 

As is well known, elastic scattering of the type a + (/3, j) m — >• a + 7)™ can 
formally be described by means of a single, one-channel LS equation of the type 

T a m,am{ z ) — ~^cmi.am( z ) Him, am ( z )Go-,am ( z )%xm,am ( z ) ■ (18) 

Though not useful for practical calculations, the very existence of an equation 
like ( |l8| ) serves as justification for use of phenomenological OP's. 

The solvability of (|l8|) depends on the properties of the plane wave matrix 
elements V°^ am (q' a , q a ; z) = (q' a \V°% iam ( z )\q a ) in the limit A a = q' a -q a -> 
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which reflects itself in a corresponding behaviour in co-ordinate space for large 
intercluster separation p a . The exact expression for the OP comes as a sum 
of two terms V°^ am = V^ am + V°^ am . The first term, so-called static 
potential, has only the trivial Coulomb-type behaviour which for a spherically 
symmetric target reads as (ep~ ( = e@ + e 7 ) 

■vistat ( i \ 47re a e/3 7 ^ tat Pa ^co e a e/3 7 

Thus it is the 'nonstatic' part V°^ Qm which is of interest. For energies 
E < 0, i.e. below the three-body threshold, its behaviour is known from 2 nd 
order perturbation theory, and has been corroborated for the exact OP in jl8| : 

Km,am(?Qi9ai-^) ~ ^a ^ r~ ^am,mn(Pa) ~ ~~ r, 4 • (20) 

Here, a static dipole polarisability of the composite particle. 

But does this result hold also above the three-body threshold? In fact, 
the two perturbative answers are in support of it, but Ref. |ll| finds an energy- 
dependent " a" while in [^0j the standard expression for a is obtained. And what 
happens to the exact nonstatic OP for E > 0? Answering requires the investi- 
gation of the zero- momentum transfer limit of am (q' al q a ; z) which is fully 
determined by the Coulomb part of the channel interaction V£ = Ylujta 
and by the three-body Coulomb resolvent G c (z) = (z — Hq — J2 U ) _1 > 

^am,am(9a! Qai z ) ~ cm jam (?a J lot j z ) 

■= (q'a\(^a m \V^Q am G c (z)Q am V^\i !am )\q a }. (21) 

Here, Q a m — 1 — IVWi) (V'aml is the projector onto all (discrete and continuum) 
target states except the incoming (= outgoing) bound state. 

The procedure goes as follows: insert the spectral representation of G c , with 
both two- and three-body scattering states; for the latter use the asymptotic 
three-charged particle wave function in Q a . Then extract the leading singularity 
in the limit A a — > 0. (Note: 2 nd -order perturbation theory is equivalent to 
replacing G c by the Coulomb channel resolvent G^ — (z — H — V^p) -1 .) 



1.2.2 Nonperturbative Result Valid for All Energies 

Let me just state the result without giving any details. For the latter I refer to 
[ pl| . Off the energy shell one finds 

Ktjal (q' a ,q a ;E + i0) Aa =° CA a +o(A a ), (22) 

with an energy- and momentum-dependent 'strength factor' C = C{E, q a ). On 
the energy shell, i.e. for E — q^/2M a + E am , q' a — q a , [E am is the bound state 
energy), Vamjam(A\ a ) is a function of the momentum transfer only. Hence, in 
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co-ordinate space one indeed ends up with a local, energy-independent potential 

V3£$(Pa) Pa =°° + ( 23 ) 

A further important point is that the exact 'strength factor' a coincides with 
the polarisability as found in perturbative approaches, 

x - \p a . D nm \ 2 f dk° a \p a -D k o a J 2 

\\Eam\ - \Ecm\\ J { } \\E am \ + /2fl a \ 

/m^-ep/mp) J dr a ip* an (r a)T a ll) am (v a ) (25) 
-DfcOm = taf)e a Pa (e 7 /m 7 - ep/mp) J dr a i}j^*{r a )r a ijj am {r a ) (26) 

being the matrix elements of the dipole operator between the incoming and all 
other (bound and continuum) target states. That is, no " renormalisation" of a 
arises from summing up all higher order terms in the perturbation expansion of 
G c , and the E- and ^-dependence of the 'off-shell strength' C has disappeared. 

Concluding I mention that a new nonrelativistic contribution to the asymp- 
totic part of the OP ~ which results from the long-range three-body 
correlations described in Sect. ( |l.l.3| ), has been derived in |p2[ . 



2 Calculational Developments 

Considerable progress has been made in recent years also in applying few-body 
theory to the practical computation of atomic and nuclear reactions. Here, I 
will concentrate on describing an attempt at re-establishing this method as a 
powerful tool for atomic physics calculations. Advances in the calculation of 
proton-deuteron scattering with realistic nuclear potentials will only briefly be 
touched upon (cf. the separate Contribution to this Conference |23|). 

2.1 Few-Body Approach to the Energetic Atomic Collision H + + H 

Scattering of protons off hydrogen has attracted the interest of physicists since 
the '30. Because direct scattering, electron exchange (both including excita- 
tion) and ionisation are coupled by unitarity they should in principle constitute 
an 'ideal' field of application of three-body theory. However, except for a few 
isolated, and not so successful, attempts little has been accomplished in this 
respect. 

Instead, most efforts have concentrated on developing, and improving, stan- 
dard models based on the SE. For instance, at low and medium energies the 
expansion into target states ('close coupling') appears to be the method of 
choice. At higher energies (from keV to a few hundreds of keV), which is the 
regime I will concentrate on, DWBA-type models or models based on the mul- 
tiple scattering expansion of Lippmann-Schwinger equation for wave functions 
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are rather successful. However, most of these approaches treat only one channel 
explicitly, i.e., intermediate-state rearrangement and/or excitation is ignored. 
One consequence is the violation of (three-body) unitarity constraints. Hence, 
though they generally provide a good description of either electron exchange 
or (in)elastic cross sections, they are frequently unable to describe the other 



reaction. For recent reviews I refer to |24|, 25 . 

Two-fragment processes of the type a + (/3, 7)™ — ► (3 + (a, j) n are conve- 
niently described by means of the effective-two-body formulation of the three- 
body theory There, the amplitudes for the various two- fragment processes 
satisfy multichannel Lippmann-Schwinger-type equations 

(XPn I G I >+ E (X(3n\G T^G \ Xa m) + ---,(28) 

with 5p a — 1 — 5/3 a . The 'form factor' | Xam) is a certain off-shell extrapolation 
of the bound state wave function | ijj am )- Already the first-order terms ('trian- 
gle amplitudes') in the expansion of V contain matrix elements of an operator 
T^, the most awkward part of which is the T-matrix describing Coulomb 
rescattering of the projectile off each of the target particles. Due to the singular 
nature of the latter, calculation of such terms is very difficult (this is particu- 
larly true if the rescattering particles have charges of opposite sign, because of 
the infinity of bound states). Indeed, the lack of applications of Faddeev-type 
approaches to atomic reactions is generally attributed to this fact |24| ]. 

Hence, the so-called Coulomb-Born Approximation (CBA) is usually made 
which consists in approximating by the Coulomb potential V^p, v = a, 0, 7. 
Recently, however, we have succeeded in calculating exactly all (on-shell) tri- 
angle amplitudes with the full (attractive and repulsive) Coulomb T-matrix, 
for all energies and scattering angles, for arbitrary bound state wave functions. 
In this way we could show that the CBA almost always fails badly for atomic 
triangle amplitudes [^7| |H, [2{|. (I mention that we have also derived (semi- 
analytical approximations for these triangle amplitudes, valid for the higher 
energies considered here, which greatly simplify practical calculations but are 



much more accurate than the CBA 129, 30 j) 



The exact triangle amplitudes have been employed in an investigation of 
direct and exchange scattering of H + off hydrogen atoms. For this purpose 
we have written Eqs. ( p7| ) in K-matrix approximation. When the resulting 
two-dimensional integral equations are transformed to the impact parameter 
representation they lead to a system of coupled algebraic equations |3l| . The K- 
matrix has been approximated by the effective potential ( p8| ) were we keep the 
th and the full 1 st order terms. For the latter, numerically calculated quantities 
the transformation to the impact parameter representation has been performed 
by quadrature. Our results provide a simultaneous, satisfactory description of 
differential cross sections for both elastic scattering and electron exchange; 



likewise the total exchange cross section is well- reproduced 1 32 
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2.2 Proton-Deuteron Scattering in the Integral Equations Approach 

The motivation behind 30 years of investigation of nucleon-deuteron (Nd) scat- 
tering rests on the hope that a detailed understanding may lead to deeper in- 
sights into the properties of 27V and 3N forces etc. Experimentally the proton- 
induced reaction (pd) is preferred; but reliable and sophisticated calculations 
have become available only recently, for energies below the deuteron breakup 
threshold |53[ Some of the underlying reasons are pointed out in J23|. Of 
course, the common practice of comparing nd calculations with pd data, though 
in praxi rather successful, is clearly unsatisfactory from principle point of view. 

We |23| have obtained the first reliable results for pd scattering above the 
deuteron breakup threshold, for the Paris potential. We used the screening 
and renormalisation approach |35j ]. Comparison with elastic scattering cross 
section and polarisation data shows a reasonably good absolute agreement; but 
a much better reproduction is achieved of the difference between experimental 
pd and nd observables where available. Further progress requires, in particular, 
excellent, low-rank, separable representations of more modern potentials. 

Acknowledgement. Most of the results presented here have been worked out in 
collaboration with G.V. Avakov, B.F. Irgaziev, A.S. Kadyrov, A. M. Mukha- 
medzhanov, and A. I. Sattarov. 
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